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We investigate the critical behaviour of charged and rotating AdS black holes in d spacetime 
dimensions, including effects from non-linear electrodynamics via the Born-Infeld action, in an 
extended phase space in which the cosmological constant is interpreted as thermodynamic pressure. 
For Reissner-Nordstrom black holes we find that the analogy with the Van der Walls liquid-gas 
system holds in any dimension greater than three, and that the critical exponents coincide with 
those of the Van der Waals system. We find that neutral slowly rotating black holes in four space- 
time dimensions also have the same qualitative behaviour. However charged and rotating black 
holes in three spacetime dimensions do not exhibit critical phenomena. For Born-Infeld black holes 
we define a new thermodynamic quantity B conjugate to the Born-Infeld parameter b that we call 
Born-Infeld vacuum polarization. We demonstrate that this quantity is required for consistency of 
both the first law of thermodynamics and the corresponding Smarr relation. 
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I. INTRODUCTION 

The study of black hole thermodynamics has come a 
long way since the seminal start by Hawking and Page 
PI, who demonstrated the existence of a certain phase 
transition in the phase space of a Schwarzchild-AdS black 
hole. Specific interest developed in the thermodynamics 
of charged black holes in asymptotically AdS spacetimes, 
in large part because they admit a gauge duality descrip- 
tion via a dual thermal field theory. This duality descrip- 
tion suggested that Reissner-Nordstrom AdS black holes 
exhibit critical behaviour suggestive of a Van der Waals 
liquid gas phase transition |2|, |3[ . 

Recently this picture has been substantively revised 
HI- By treating the cosmological constant as a thermo- 
dynamic pressure and its conjugate quantity as a ther- 
modynamic volume, the analogy between 4-dimensional 
Reissner-Nordstrom AdS black holes and the Van der 
Waals liquid-gas system can be completed, with the crit- 
ical exponents coinciding with those of the Van der Waals 
system and predicted by the mean field theory. While 
previous studies considered the conventional phase space 
of a black hole to consist only of entropy, temperature, 
charge and potential, investigations of the critical be- 
haviour in the backdrop of an extended phase-space (in- 
cluding pressure and volume) seems more meaningful for 
the following reasons. 

• With the extended phase space, the Smarr relation 
is satisfied in addition to the first law of thermody- 
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namics, from which it is derived from scaling argu- 
ments. This makes it evident that the phase space 
that should be considered is not the conventional 
one @- 

• The resulting equation of state can be used for com- 
parison with real world thermodynamic systems [3- 
|8(. For example, the variation of these parameters 
makes it possible to identify the mass of the black 
hole as enthalpy rather than internal energy. As 
a consequence, this may be useful in interpreting 
the kind of critical behaviour (if any) to the known 
thermodynamic systems (if any). 

• Thermodynamic volume has been studied for a 
wide variety of black holes and is conjectured to 
satisfy the reverse isoperimetric inequality Q. 

• Variation of parameters in the extended phase 
space turns out to be useful for yet another reason: 
one can consider more fundamental theories that 
admit the variation of physical constants [H, Q . 

As noted above, a recent investigation of 4-dimensional 
charged (Reissner-Nordstrom) AdS black hole thermody- 
namics in an extended phase space — where the cosmolog- 
ical constant is treated as a dynamical pressure and the 
corresponding conjugate quantity as volume — indicated 
that the analogy with a Van der Waals fluid becomes very 
precise and complete One can compare appropri- 
ately analogous physical quantities, the phase diagrams 
become extremely similar, and the critical behaviour at 
the point of second order transition identical. 

In this paper we elaborate on this study by considering 
more general classes of charged and rotating black holes: 
specifically, the d-dimensional analogues of the Reissner- 
Nordstrom- AdS black hole, its nonlinear generalization 



2 



described by Born-Infeld theory, and the Kerr-Newman 
solution. We find that for d > 3 the analogy with the Van 
der Waals fluid holds, with the parameters appropriately 
generalized to depend on the dimension. However we 
show that for d = 3 neither the charged nor rotating BTZ 
black holes exhibit critical behaviour. While we find the 
critical exponents in all cases to be the same as for the 
4-dimcnsional case, we uncover interesting new thermo- 
dynamic properties in the Born-Infeld case. Specifically, 
we find a term conjugate to the Born-Infeld parameter 
that we interpret as the Born-Infeld polarization of the 
vacuum. We correspondingly enlarge the phase space, 
allowing the Born-Infeld parameter (and its conjugate) 
to vary. This leads to a complete Smarr relation satisfied 
by all thermodynamic quantities. 



II. PHASE TRANSITION OF CHARGED ADS 
BLACK HOLES IN HIGHER DIMENSIONS 

A. Charged AdS black holes 

The solution for a spherical charged AdS black hole in 
d > 3 spacetime dimensions reads 



ds z 



-Vdt 1 



dr 2 



r 2 dtf d _ 2 , 



F = dA. 



A 



1 d-2 
2(d-3) 



r d— 3 



dt. (2.1) 



Here, dSl 2 d stands for the standard element on S and 
function V is given by 



V = 1- 



The outline of our paper is as follows. We first look 
at higher dimensional Rcissner-Nordstrom black holes to 
obtain the dependence of the various phase space quan- 
tities and critical exponents on the spacetime dimension- 
ality d. We find that the critical exponents do not de- 
pend on the dimension provided d > 4. We then con- 
sider the Kerr-Newman solution and find in the limit of 
fixed small angular momenta J that the the same sit- 
uation holds qualitatively, with the same exact critical 
exponents. Turning to the d = 3 case, we then consider 
charged BTZ black holes. In this case there is no critical 
phenomena, a situation that persists even if we include 
rotation. We then investigate non-linear electrodynamics 
via the Born-Infeld action in four spacetime dimensions. 
We encounter a rich structure when we investigate the 
Born-Infeld black hole case. We can have a singularity 
cloaked by one or two horizons or have a naked singular- 
ity depending on the value of the parameters. We find 
that in order to satisfy the Smarr relation we need to 
consider the Born-Infeld parameter as a thermodynamic 
phase space variable; in so doing, we introduce its con- 
jugate quantity, polarization, and calculate it. We com- 
pute the appropriate Smarr relation and then study the 
critical behaviour. We close with a concluding section 
summarizing our results. Appendix A reviews the defi- 
nition of critical exponents and their values for the Van 
der Waals fluid and the Reissner-N6rdstrom-AdS black 
hole. In appendix B we present an heuristic derivation 
of the Van der Waals equation in higher dimensions. We 
gather in appendix C various hypergeometric identities 
used in the main text. In what follows we always assume 
Q > without loss of generality. 
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(2.2) 



Parameters m and q are related to the ADM mass M (in 
our set up associated with the enthalpy of the system as 
we shall see) and the black hole charge Q as 0] 



M 



d-2 
16tt 



■Wrf_ 2 TO, 



y/2(d-2)(d-3) 

Q = 5 q , 

with Ud being the volume of the unit d-sphcrc, 



2-7T 2 



(2.3) 



(2.4) 



The metric and the gauge field (|2.1j) solve the Einstein- 
Maxwell equations following from the bulk action 



(d- l)(d- 2)> 




In our considerations we interpret the cosmological con- 
stant A = — (£L_U^L_E2 as a thermodynamic pressure P, 



P 



1 

8^" 



-A 



{d- l)(d- 2) 
16ttZ 2 



(2.6) 



The corresponding conjugate quantity, the thermody- 
namic volume, is given by Q 



V 



d-l 

The black hole temperature reads 



(2.7) 



Note added. We note here recent work concerned 
with angular momentum that has some overlap with our 
paper |10| . 



T 



V'(r+) _d-3 
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2(d-3) d _ 3 12 



(2.8) 



with 7+ being the position of black hole horizon, deter- 
mined from V{r + ) = 0. The black hole entropy S and 
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the electric potential $ (measured at infinity with respect 
to the horizon) are 



g — Ad-2 A _ , , . „d-2 



d-2 q 



(2.9) 
(2-10) 



y 2(d - 3) r^ 3 ' 
All these quantities satisfy the following Smarr formula: 



M=^^TS + ^Q--^—VP, 
a — 3 a — 3 



(2.11) 



as well as the (extended phase-space) 1st law of black 
hole thermodynamics 



dM = TdS + $dQ + VdP . 



(2.12) 



These two are related by a dimensional (scaling) argu- 
ment; see, e.g., @. 

In the canonical ensemble (fixed charge Q) a first order 
phase transition in the (T, <5)-plane, reminiscent of the 
liquid-gas phase transition of the Van der Waals fluid 
was noted quite some time ago in the context of a fixed 
cosmological constant 0, H| . The properties of the corre- 
sponding critical point have been elaborated, for exam- 
ple, in [11(. More recently, however, it has been pointed 
out [|| that at least in 4D the corresponding critical be- 
haviour has a more natural interpretation in the extended 
phase space (where the cosmological constant is treated 
as a thermodynamic pressure), giving the phase transi- 
tion in the (P, T)-plane, making the similarity with the 
Van der Waals fluid complete. In what follows we shall 
study this situation in the higher-dimensional case. 



The behaviour of G for dimensions d = 4 and d = 10 
is depicted in Figs. [TJ and [2J respectively; we see from 
these figures the characteristic lst-order phase transition 
behaviour. 

To make contact with the Van der Waals fluid in d 
dimensions (see App. B), note that in any dimension d 
we have lp~ 2 = G c ih/c 3 , where pressure has units of en- 
ergy per volume. We thus identify the following relations 
between our 'geometric quantities' P and T and the phys- 
ical pressure and temperature: 



[Press] = 

Therefore 

Press 



(2.16) 
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(2.17) 
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Comparing with the Van der Waals equation, (|B9|) . we 
conclude that we should identify the specific volume v of 
the fluid with the horizon radius of the black hole as 



Ar+l 



d-2 



d-2 



In geometric units we have 



d-2 



(2.18) 



(2.19) 



B. Critical behaviour 

The critical behaviour of a system is captured by its 
partition function. Namely the thermodynamic poten- 
tial, associated with the the Euclidean action, calculated 
at fixed Q, fixed P and fixed T is the Gibbs free en- 
ergy. Using the counterterm method (canceling the AdS 
divergences) [H, [l3[ and the following boundary term 
(characteristic for the canonical ensemble [HI): 



1 

-in 



h = — / d?- 1 xVhF ab n a A b , (2.13) 



the Gibbs free energy G = M — TS reads @ 



G = G(P,T) 



u d-2 

16tt 



VotxPt 



d-i 



{2d-b)q' 



+ 



(d-l)(d-2) r 



d-3 



(2.14) 

Here, r+ is understood to be a function of the black hole 
temperature T and pressure P via equations (|2.6p and 
(I2.8[) : upon employing <\2.7\ we obtain 



P 



T(d-2) (d-3)(d-2) q 2 (d-3)(d-2) 



4r 4 



167T7 + 



16nr 



2(d-2) 
+ 



(2.15) 



and the equation of state reads 



P 



T 



(d-3) 



q 2 (d-3) 



n{d-2)v 2 ^ v 2{d-2) K 2d-5 



(2.20) 



The associated P—V diagrarrQ for a 4D black hole is dis- 
played in Fig. [3] Obviously, for T < T c there is a small- 
large black hole phase transition in the system. This 
qualitative behaviour persists in higher dimensions — see 
Fig-Hlfor an illustration in d = 10. 

Critical points occur at points of inflection in the P — V 
diagram, where 



dP _ d 2 P 
dv dv 2 



0, 



(2.21) 



1 Throughout our paper we display P — v diagrams, with v being 
the specific volume of the corresponding fluid, rather than P — 



V diagrams, with V being the thermodynamic volume, 
dimensions the two are related via 
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FIG. 1. Gibbs free energy of charged AdS black hole 

in d = 4. Gibbs free energy G of 4D RN-AdS black hole is 
depicted as function of temperature T for fixed q = 1 and 
various pressures P/P c = 1.6,1,0.6 and 0.2. The behaviour 
for d > 4 is qualitatively similar, see next figure. 
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FIG. 2. Gibbs free energy in d — 10. Gibbs free energy G 
of 10D RN-AdS black hole is depicted as a function of temper- 
ature T for fixed q = 1 and pressures P/P c = 1.6, 1, 0.6, 0.2. 




FIG. 3. P — V diagram of charged AdS black hole in 

d = 4. The temperature of isotherms decreases from top to 
bottom. The two upper dashed lines correspond to the "ideal 
gas" one-phase behaviour for T > T c , the critical isotherm 
T = T c is denoted by the thick solid line, lower (red) solid 
lines correspond to two- phase state occurring for T <T C . We 
have set q = 1. The behaviour for d > 4 is qualitatively 
similar - see figure 2] 




FIG. 4. P — V diagram in d = 10. The description coincides 
with that of the previous figure. 



yielding 



v c = - [q 2 {d-2)(2d-5] 



Pr 



l/[2(d-3)] 



(rf-3) 2 
itnv c (2d — 5) 

(cZ-3) 2 
167rK 2 w c 2 



(2.22) 



This gives the following universal (independent of q) ratio 
P c v c 2d - 5 



T, 



4d 



(2.23) 



Note that only in d = 4 do we recover the ratio p c = 3/8, 
characteristic for a Van dcr Waals fluid in any number of 
spacetime dimensions. 
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0.004 - 



Critical Point 



izing the behaviour of physical quantities in the vicinity 
of the critical point found in the previous section. The 
definition of these exponents is for convenience gathered 
in appendix A. 

To calculate the critical exponent a we consider the 
entropy S, ([2T9"]| . as a function of T and V. Using (f2~7) 
we have 
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FIG. 5. Coexistence line of charged AdS black hole in 

d = 4. This figure displays the coexistence line of the small- 
large black hole phase transition of a charged AdS black hole 
for d = 4 system in the (P, T)-plane. The critical point is 
highlighted by a small circle at the end of the coexistence line. 
The behaviour for d > 4 is qualitatively similar, as shown in 
figure [6] 
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FIG. 6. Coexistence line in higher dimensions. The 

coexistence line is displayed for d = 4, 5, 6 and 10. 



The coexistence line in the (P, T)-plane is displayed in 

Figs, na 



C. Critical exponents 

Following the approach taken previously for d = 4 Q , 
we shall now calculate the critical exponents character- 



s' = S(T, V) = 



'd-2 



(d-l)V 



(2.24) 



Since this is independent of T, we have Cy = and hence 
a — 0. Defining 



P= P, 



v 

Vr. 



T 



(2.25) 



the equation of state (|2.20j) translates into the following 
'law of corresponding states': 



V 



4(d- 2)r 



d-2 



1 



(2d-b)v (d-3)v 2 (d- 3) (2d- 5)is 2d ~ 4 ' 

(2.26) 

Although this law has explicit dimensional dependence, 
we shall show that this does not affect the behaviour of 
the critical exponents. 

Slightly more generally, let us assume (as will be the 
case for all the examples studied in this paper) that the 
law of corresponding states takes the form 

1 r 

Pc V 



P 



(2.27) 



where p c stand for the critical ratio. Expanding this 
equation near the critical point 



T = t+1, v=(u + l) 



1/9 



(2.28) 



where q > 0, and using the fact that from the definition 
of the critical point we have 



Pc 

and so obtain 



1 



Pc/'(1) = 1, Pc/"(l) = -2, (2.29) 



p=l + At- Btu - Cuj 3 



where 



A 



1 

Pc 



D 



1 

qp c 



c 



0(icj 2 ,ij 4 ), (2.30) 
/(3) (1) - 



6 



(2.31) 

Let us further assume that C > (again this is will be 
the case of all the examples). Differentiating the series 
for a fixed t < we get 

dP = -P c (Bt + 'iCuj 2 )duj . (2.32) 

Employing Maxwell's equal area law, see, e.g., [H, while 
denoting uj s and uji the 'volume' of small and large black 
holes, we get the following two equations: 



p = 1 + At - Btuii - Cujf = 1 







codP. 



At - Btuj s - CujI 

(2.33) 
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The unique non-trivial solution is 



The U(l) potential reads 



uj s = —u>i = 



-Bt 
~C~ 



yielding 

rj = V c [u>i —u) a ) = 2V c ui oc \f^i 



(2.34) 



(2.35) 



To calculate the exponent 7, we use again (|2.32j) . to get 

1 1 



Kt 



l_dV 
VdP 



T ^ P r Bt 



7 



1 



(2.36) 



Finally, the 'shape of the critical isotherm' t = is given 
by U23DJ), i.e., 



p - 1 = -Cuj 3 



6 = 3. 



(2.37) 



To conclude, provided that the equation of state takes 
the form (|2.27[) . such that the coefficient C, given by 
(|2.3ip is non-trivial and positive, we recover the (mean 
field theory) critical exponents 



7 = 1, 5 = 3. 



(2.38) 



In particular, for charged AdS black holes in any di- 
mension d, the law of corresponding states, (|2.26[) . takes 
the form (|2.27|) . Taking q = d — 1 (in which case 
ui = X- — 1), we obtain the expansion (|2.30|) with 



A = 



Ad- 8 
2d - 5 



B 



4d-8 



(2d-5)(d- 1) ' 



C 



2d- A 



3(d- l) 3 
(2.39) 

and so the discussion above applies. We conclude 
that the thermodynamic exponents associated with the 
charged AdS black hole in any dimension d > 3 coincide 
with those of the Van der Waals fluid. We now consider 
the effects of rotation on these results in the following 
section. 



III. ROTATING BLACK HOLES 

A. Thermodynamics 

The charged AdS rotating black hole solution is given 
by the Kerr-Ncwman-AdS metric 



= 4 



dt 



S sin 2 9 



a sin 9 
— ^ — dip 



adt 



-dip 



d = 4, where 



p 2 = r 2 + a 2 cos 2 £ 



,.2 ' 



I 2 ' 



5=1 



(3.1) 



■ cos 2 9 . 



A = (r + a ) ( 1 + — ) - 2mr 



(3.2) 



or / , asm u , 
A = -^r [ tit — dip 



(3.3) 



The thermodynamics of charged rotating black holes 
treating the cosmological constant as independent vari- 
able was performed in @, BHH- In particular, one has 
the following thermodynamic quantities: 



s _ ir(r 2 + + a 2 ) 



T = 



-1 1 a 2 1 o r + a 2 +q 2 



A-n(r\ + a 2 ) 



rl + a 2 



r\ + a 2 



(3.4) 



The mass M, the charge Q, and the angular momentum 
J are related to parameters to, q, and a as follows: 



TO _ o T aTO 



(3.5) 



Using the counterterm method [12|, Il3j , th e action for the 
canonical ensemble was calculated in [ljl [l6| and reads 



/3 ( r\ a 2 + q 2 2q 2 r + 



1 4 V + 3 2 ' r + E ' S(r 2 + a 2 ) 
The corresponding Gibbs free energy reads 



G = G(P,T,J,Q) = - + flJ, 



where [if 



/ 2 



The thermodynamic volume is [8|, Q 

2ir {r\ + a 2 ){2r 2 + l 2 + a 2 l 2 - r\a 2 ) + l 2 q 2 a 2 



V 



l 2 ~ 2 r + 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



with the pressure P = -^j? still given by (|2.6[) . The 
equation of state is written as 



T 



1 -J- ° 2 4- t r + _ ° 2 +9 2 



47r(r^ + a 2 ) 



(3.10) 



where we express the parameters a and q in terms of the 
physical quantities J and Q, (|3.5|) . This gives 



P = 



r 1 g 2 



2r + 87^ %irr\ 

3(4r| + 87rTr5 - 87rT?' 3 Q 2 + r^Q 2 - 2Q 4 ) 2 
+ 87rr6 (r^.+ 2^V| + 2Q 2 ) 2 J 

+0(J 4 ). (3.11) 

as an expansion in powers of the angular momentum J. 
In what follows we expand all quantities to 0(J 2 ), ne- 
glecting all terms higher order in J. We further introduce 



7 




/ 



2 4 /'6 8 10 12 14 16 18 20 



FIG. 7. P — V diagram for the rotating black hole. The 

isotherms in P — V plane are depicted for a rotating black 
hole with J = 1. 
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FIG. 8. Gibbs free energy for the rotating black hole. 

Isobars of the Gibbs free energy for a rotating black hole with 
J = 1 are depicted for P/P c = 1.6, 1, 0.6, 0.2. 



the following quantity (which we shall see corresponds to 
the specific volume of the associated Van der Waals fluid) 



1/3 



2r. 



12(87rr-4p 



3rl 



4tt J '~ r ' r+(3r 2 + %irr\P ■ 

and finally obtain the equation of state 



3Q 2 ) 



2\2 



(3.12) 



T 1 
P = =• 

v 2ttv 2 



2Q 2 48J 2 



7TV 



96Q 2 (24Q 2 + 5v 2 + 6itTv 3 )J 2 
+ v 2 + nTv 3 ) 2 



(3.13) 



In what follows we limit ourselves to the case of zero 
charge Q = 0, for which (|3.13[) becomes Q 



P 



T 

v 



1 



48 J 2 



2ttv 2 



(3.14) 



In the same approximation we get the following expres- 
sion for the Gibbs potential: 



v irPv* 
G= 8~ ~Y2~ 



20J 2 



(3.15) 



Upon comparing (f3~14|) to (f2T20|) and (|3~T5| to (|2~14)) we 

note the same qualitative behaviour as in the charged 
case, with fixed J playing the role of fixed Q. However 
the quantitative dependence of P on v is different. We 
depict the P — V diagram in Fig. [7] and the Gibbs poten- 
tial in Fig. [SJ where the qualitative behaviour is seen to 
be similar to the non-rotating charged case. The P — T 
diagram is illustrated in Fig. |H1 
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FIG. 9. Coexistence line for the rotating black hole. 

The coexistence line of the phase transition between small and 
large slowly rotating black holes. The figure is qualitatively 
similar to charged AdS black hole case. 



Critical behaviour 



Using (|2.2ip the critical point occurs at 



Pr 



2 x 90 1/4 \/j, 
90 3 / 4 1 
2257 VJ ' 
1 1 



12v/907r J 



(3.16) 
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and these relations in turn yield the universal critical 
ratio 



P 



Pc = 



5 

12 



(3.17) 



that differs from the Van der Waals fluid. Upon inclusion 
of higher powers of J we expect this to become J depen- 
dent. We shall see a similar phenomenon take place in 
the Born-Infeld case. 

Defining the quantities p, v and r as previously, (|2.25l) . 
we get the following law of corresponding states: 



12 r 3 1 



1 1 



P 5 v 2 v 2 10 1/ 6 



(3.18) 



which is again universal and of the form (|2 . 27"[) . Expand- 
ing around the critical point using the variables r and u> 
with f/ = 3, we get 



p= H i tw lj 3 + 0(tuj 2 ,u> A ). 

P 5 5 27 V ' ; 



(3.19) 



It is now obvious from this expansion and the discussion 
in Sec. II. C that the critical exponents /3, 8 and 7 do not 
change, i.e., we have 



Moreover, we have 



7=1, 6 = 3. 



(3.20) 



S = S(T, V) = 



ir{r 2 



TTV 



0(J 4 ). (3.21) 



Hence to quadratic order in J we have 



a, 



dT 



0, 



(3.22) 



which implies that the critical exponent a = 0. So we 
conclude that (at least for slowly rotating black holes) 
the critical point belongs to the same universality class 
as the one of the Van der Waals fluid. 

One might wonder if this continues to hold also for 
d = 3. We show in the next section that this is not the 
case, and the BTZ black hole is rather different from its 
charged/rotating counterparts in higher dimensions. 



with 



Q 2 /r\ r 2 



(4.2) 



With this choice we have a solution of the Einstein- 
Maxwell system with the cosmological constant A = 
— l/l 2 = — 8irP. The temperature is 



T = 



V(r+) 



4tt 



2ttZ 2 



(4.3) 



which translates to the following equation of state: 



P = 



T 



2ttv 2 



v = 4r 'a 



(4.4) 



It is obvious that such an equation does not admit any 
inflection point and hence the charged BTZ black hole 
does not exhibit any critical behaviour. This is quite 
analogous to the toroidal (planar) case observed in four 
dimensions Q. 

One might wonder whether this remains true when one 
includes rotation. The corresponding metric was con- 
structed in [l8| and is rather complicated. Let us limit 
ourselves on considering the uncharged rotating BTZ 
black hole. The rotating BTZ black hole is given by [ItJ 



ds 2 



Vdt 2 



J 



dt + dip 



dr 2 

+ ^2 ■ 

The corresponding equation of state reads 



V = -M 



P 



T 8J 2 



4r 4 



(4.5) 



(4.6) 



Similar to the charged case it does not admit critical 
behaviour. We expect this to remain true also in the 
general rotating case with nontrivial charge (l8j . 

To conclude, for BTZ black holes, contrary to the 
higher-dimensional case, we do not observe any critical 
behaviour. 



V. INFLUENCE OF NONLINEAR 
ELECTRODYNAMICS 



IV. BTZ BLACK HOLES 



Born— Infeld-AdS solution 



In this section we investigate the case of a charged BTZ 
black hole. This case needs to be considered separately 
due to the behaviour of the electromagnetic field in three 
spacetime dimensions. The metric and the gauge field 
are given by [13 



ds 2 = -Vdt 2 4 
F = dA, A 



dr 2 

— + r 2 d V 2 . 



-Qln(-)dt 



(4.1) 



The Einstein-Born-Infeld theory in AdS is described 
by the following bulk action [l9| 



Iem = -— [ \[ ZI g \R + C B i + 4 
16tt J m \ V 



£ B / = 4& 2 (1- a/1 



2F 



F=- A F ab F ab (5.1) 



where we restrict ourselves to d = 4. The parameter 
b represents the maximal electromagnetic field strength. 



FIG. 10. BI electric field strength. The field strength E 
is depicted as a function of r for various values 6 and fixed 
Q — 1. The value of 6 decreases from top to bottom. The 
lower three red solid lines correspond to b = 0.2, b = 1/2, and 
6=1, respectively. The field reaches the finite value in the 
origin in these cases. The top (black) dashed line corresponds 
to the limiting Maxwell case (6 — > oo). 



FIG. 11. Types of Born-Infeld-AdS black holes. Two 

types of possible BI-AdS black hole solutions are displayed for 
6 = 2. The marginal case M = Mm is highlighted by a thick 
solid line. The (red) solid lines below it represent the S-type 
of M > M m - The upper dashed lines correspond to various 
cases of RN-type. The naked singularity, extremal BH, and 
two horizon solution of this type are displayed from top to 
bottom. The parameters were chosen Q = 1 and 2 = 1. 



With motivations from string theory, e.g. [2fJ, b can also 
be related to string tension, b = tt—t ■ 

In Schwarzschild-like coordinates the metric and the 
electroma gne tic field of the spherically symmetric solu- 
tion read [2l|-|i 



dv 

ds 2 = -Vdt 2 + + r 2 dnl , 

v 

F = Edt A dr , E 



Q 



^r 4 + Q 2 /b 2 



(5.2) 
(5.3) 



Here, df2| stands for the standard element on S 2 and 
function V is given by 



V=l- 



2M 

r 
2M 



4Q^ 
3r 2 



r 2 2tf_ 
l 2+ ~ 

r 2 2b 2 r 2 

F ( l l5 - 

2^1 I 7) oi 7' ~~ 



— -r 2 )dr 



4' 2' 4' bh 



(5.4) 



with 2-Fi being the hypergeometric function. Note that 
whereas the hypergeometric series expression is conver- 
gent only for \z\ < 1, i.e., for r > tq = y/Q/b, one may 
use the integral representation, (|C6j) (discussed further 
in the appendix) for any r > 0. 

The parameter M represents the ADM mass and pa- 
rameter Q the asymptotic charge of the solution. The 
electric field strength is depicted in Fig. [TO] The struc- 
ture of the horizons is discussed in the next subsection. 



B. Two types of black hole solutions 



The solution ()5.2[) possesses a singularity at r = 0. 
This singularity is cloaked by one, or two horizons, or 
describes a naked singularity, depending on the value of 
parameters. There are two types of black hole solutions 
as sketched in Fig. [TTJ This can be seen from the expan- 
sion of V around r = 0. Using the formula (|C9[) derived 
in App. C, we find (HH 



V = 1 - 



2(M-M m ) 



-2bQ + 0(r), (5.5) 



where 



M m = 



6 



(5.6) 



is the 'marginal mass '. Depending on the value of M, b, 
I, and Q we have the following cases: 

For M > M m we have the 'Schwarzschild-like' (S) type. 
This type of black hole is characterized by the existence 
of spacelike singularity and possesses one horizon. The 
characteristic behaviour of V is displayed in Fig. [TT] by 
(red) solid lines. The 'Reissner-Nordstrom' (RN) type 
is characterized by M < M m . Similar to the Reissner- 
Nordstrom solution, this may have zero, one, or two hori- 
zons, see Fig. [TT] and Fig. [12] The case M = M m is the 
'marginal' case, for which function V m — V(M m ) ap- 
proaches the (finite) value V m (0) = 1 — 2bQ. When this 
is positive, i.e., for bQ < i , the RN phase describes a 
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b=0.01 








b=j/g_— 
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b=4 

















FIG. 12. Marginal case. The sequence of marginal metric 
functions Vm, is depicted for various values of b and fixed Q — 
1,2 = 1. As b decreases the marginal line moves upwards, 
shifting the V of the RN-type to more positive values. For 
bQ < 1/2 the corresponding V is necessarily positive and 
hence the RN-type describes a naked singularity. 
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FIG. 13. Black hole horizon radia. The horizon radii are 
depicted for various types as a function of M/M m . Whereas 
for bQ < 1/2 only S-type black hole exists (lower red curve), 
for b > 1/2 the S-type smoothly joins the RN-type (black) 
and continues till extremal mass M ex /M m . The lower black 
dotted line corresponds to the position of inner horizon of 
RN-type. The units were chosen Q = 1 , 2 = 1. 



naked singularity and the only possible black hole solu- 
tion is the S-type for M > M m . On the other hand when 



bQ> 



(5.7) 



the RN-type describes a black hole for M ex < M < M m , 
with M ex being the mass of the extremal RN-type black 
hole, determined from V(r = r ex ) = 0, with r ex given by 



1 



2b 2 



- 2b^/rj x b 2 + Q 2 = . 



(5.8) 



We gather the possibilities for BI-AdS black holes in the 
following table, the corresponding dependence of r + on 
M/M m is depicted in Fig. [H 



BI-AdS black holes 


Condition(s) 


Type 


# horizons 


M > M m 


S 


one 


M ex < M < M m 


RN 


two 


bQ>\ 






M = M ex < M m 


RN 


one 


bQ>\ 


(extremal) 





(5.9) 



In all other cases we have a naked singularity. Let us 
stress that M m depends on b and Q, whereas M ex = 
M ex (b, Q, I). Therefore, for example, M ex < M m imposes 
a nontrivial restriction on the parameters of the solution. 
We also emphasize that the possibility of S-type BI-AdS 
black holes is sometimes completely ignored in a consid- 
eration of BI-AdS black hole thermodynamics [25|. We 



shall see that it is the S-type black holes for which the 
most interesting behaviour occurs. 



C. First law, Smarr formula, and vacuum 
polarization 

The Smarr relation follows from the first law of black 
hole thermodynamics and a scaling dimensional argu- 
ment @. [See also ^ for a generalization to the Love- 
lock gravities.] Beginning with Euler's theorem (for sim- 
plicity formulated only for two variables): given function 



/(x, y) such that f(a p x, a q y) 



rf(x, V) = px 



'21 

dx 



a r /(x, y), it follows that 

df 
dy, 



qy 



(5.10) 



To obtain the proper scaling relations, for the Born 
Infold case we must consider the mass of the black hole 
M to be a function of entropy, pressure, angular mo- 
mentum, charge, and Born-Infeld parameter b, M = 
M(S, P, J, Q, b). Performing the dimensional analysis, we 
find [M] = L , [S] = L 2 , [P] = L~ 2 , [J] = L 2 , [Q] = 
L , [b] = [E] = L _1 . Since & is a dimensionfull parame- 
ter, the corresponding term will inevitably appear in the 
Smarr formula. It is also natural to include its variation 
in the first law. Hence we have 



-2P(™ 
\dP 

fdM 



(5.11) 
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Defining further a quantity conjugate to b, 
B 



fdM\ 



the first law takes the forriQ 

dM = TdS + VdP + fid J + $dQ + Bdb . 



(5.12) 



(5.13) 



Using this and Eq. (|5.11|) . we find the following general- 
ized Smarr formula for stationary Born-Infeld-AdS black 
holes: 



M = 2{TS -VP + QJ) + <£Q-Bb. 



(5.14) 



Let us emphasize that even in the case when the cosmo- 
logical constant A and the parameter b are not varied in 
the first law, i.e., we have dM = TdS + QdJ + the 
Smarr relation (|5.14[) is valid and the VP and bB terms 
therein are necessary for it to hold. This resolves the 
problem of inconsistency of the first law and the corre- 
sponding Smarr relation raised by Rasheed [13] ; see also 

Hi El. 

In particular, for our (static) black hole solution (|5.2p 
we have the following thermodynamic quantities: The 
black hole temperature and the corresponding entropy 
are 



T = 



1 



S 



47rr+ 
A 

7' 



p 

A = Airr 2 , 



-2b 2 r 2 




(5.15) 



(5.16) 



The thermodynamic volume V and the corresponding 
pressure P are 



p = -ir A =irj2- ( 5 - 17 ) 



The electric potential $, measured at infinity with re- 
spect to the horizon, is 



Qdx 



Q 



-.Fx 



1 15 

4' 2' 4' 



Q 2 



. ^1 + Q 2 /b 2 r+ i -V4'2'4' b 2 r%, 

(5.18) 

Using the first law (|5.13p and the formula (|C13j) . we get 
the following expression for the quantity B: 




(5.19) 



2 In fact, this form of the first law remains to be proved. In or- 
der to do so, one could, for example, employ the Hamiltonian 
perturbation theory techniques, similar to l26t . 



Alternatively, one could calculate B from the Smarr for- 
mula (|5.14|) . getting 



n 2 h a A u Q2 U Q2 w ( l 1 5 Q2 



(5.20) 

The two expressions of course coincide and are related 
by (|C14|) . We shall return to this quantity later in this 
section. 



D. Critical behaviour 

The thermodynamic behaviour of Born-Infeld-AdS 
black holes and the corresponding phase transitions were 
studied in grandcanonical ensemble in I24I and in canoni- 
cal ensemble in 



21|. The critical exponents were recently 



calculated in In this subsection we discuss the 

thermodynamics in canonical ensemble in the extended 
phase space. We also correct the values of the critical 
exponents and some of the previous results appearing in 
the literature. In particular, we limit ourselves to the 
case when b is fixed and consider P — V extended phase 
space. Our aim is to study the influence of the nonlin- 
earity of the electromagnetic field on the existence of the 
critical point and its corresponding behaviour. 



Equation of state 



For a fixed charge Q, Eq. (|5.15|) translates into the 
equation of state for a Born-Infeld AdS black hole, P = 
P(V,T), 



P 



T 

277 




8irr 2 + 



(5.21) 



Similar to the Maxwell 4D case, discussed in Sec. II, it 
is the horizon diameter 2r+ , rather than the thermody- 
namic volume V, that should be associated with the spe- 
cific volume of the corresponding fluid; v = 2r + lj 
ing l P = 1, Eq. (|5.2ip gives 



l 2 p . Tak- 



P 



T 

v 



1 



2ttv 2 



-in 



1 



1 



16Q2 



(5.22) 



The corresponding P — V diagrams are displayed in 
Figs. [HHT7] for various values of b. 

The behaviour of the isotherms depends crucially on 
how 'deep' we are in the Born-Infeld regime. For bQ > 
1/2 we have the 'Maxwcllian regime' and the behaviour 
is qualitatively similar to that of a Rcisner-Nordstrom- 
AdS black hole. Namely, for Q ^ and for T < T c there 
is an inflection point and the behaviour is reminiscent of 
the Van der Waals g shown in Fig. [T4l 

However, for b < 1/2 we have much 'stranger be- 
haviour', with one, two, or zero inflection points, as 
shown in Figs. I15H17I . The critical point is obtained 
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FIG. 14. P-V diagram of BI-AdS black hole for 6 > b 2 . 

The temperature of isotherms decreases from top to bottom. 
The two upper dashed lines correspond to the "ideal gas" one- 
phase behaviour for T > T c , the critical isotherm T = T c is 
denoted by the thick solid line, lower solid lines correspond to 
temperatures smaller than the critical temperature. The fig- 
ure is very much reminiscent of the corresponding behaviour 
of RN-AdS black hole. For b > 63 ~ 1.72846/Q the critical 
point occurs for the RN-type black hole whereas for b < 63 
for the S-type. In both cases the behaviour remains quali- 
tatively the same. We have set Q = 1 and 6=1 for which 
T c « 0.04362. 




FIG. 15. P-V diagram of BI-AdS BH for b G (61,62). 
We have now two critical points, one at positive pressure, the 
other at negative pressure. The upper one occurs at bigger 
horizon radius, bigger mass, and higher temperature. Both 
are associated with the S-type of black hole solutions. We 
have set Q = 1 and 6 = 0.45 for which T C1 « 0.04517 and 
T C2 ~ 0.026885. As we decrease 6 further, the critical points 
start to move towards each other. At 6 = 61 ~ 0.4237/Q both 
have positive pressure, see Fig. 1161 Decreasing 6 even further, 
the two critical points eventually merge and disappear. That 
is for 6 < 60 ~ 0.35355/Q there are no critical points anymore, 
Fig El 



from Eq. (|2.21|) . which leads to the cubic equation 



-px 



0. 



where 



36 2 
32Q 2 



In order that v r = 



256Q 4 



16Q 2 
b 2 



(5.23) 



-1/2 



nicy 

6 2 



an additional constraint 



\x\ < 



4Q 



(5.24) 

be positive, we require 
(5.25) 



To write an analytic expression for v c in terms of 
trigonometric and/or hyperbolic functions let us analyze 
the cubic. The cubic of the form (|5.23[) . with p < and 
real q, admits three or one real roots. Three real roots 
occur when 4p 3 + 27q 2 < 0, i.e., when 



b> bn = 



0.35355 



V8Q Q 
in which case the roots are given by 

Xk = 2 4 — — cos — arccos —\ — 



(5.26) 



2p\ p 



2irk\ 

0,1,2. 



(5.27) 



0.005 




-0.005 



FIG. 16. P-V diagram of BI-AdS BH for 6 G (60,61). 
We have now two critical points at positive pressure. The 
upper one has higher radius, temperature, and mass. We 
have set Q = 1 and 6 = 0.4 for which T C1 w 0.045890 and 
T c , w 0.040492. 
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FIG. 17. P-V diagram of BI-AdS BH for 6 < bo. There 
are no longer any critical points below bo ~ 0.35355/Q- We 
have set Q = 1 and b = 0.2. 



It is easy to verify that the additional condition (|5.25j) is 
never satisfied for x 2 , always satisfied for x\, and satisfied 
for xq provided that 

-^-hStSt-i. (5.28) 

One real root occurs for b < bo , in which case the solution 
is given by 



x 3 = -2 



^coshfiarccoshf-^ 



2p\ p 



(5.29) 



which, however, does not satisfy (|5.25[) . 

To summarize, the critical point can occur only for b > 
&o- Between bo < b < 62 we have two critical points, given 
by (|5.27p . with x = £0,1, where x\ corresponds to a larger 
v c . Above 62 there is one critical point, given by x = x\. 
Using a series expansion, it is straightforward to show 
that for b — > 00, v c determined from x\ approaches the 
'specific volume' of the Reissner-N6rdstrom-AdS black 
hole. That is why we call the branch determined from x\ 
the Maxwell (M-branch), and the one determined from 
xq (which is deeply in the Born— Infeld regime) the BI- 
branch. Once the critical radius v c has been determined, 



1 16Q 2 \ J 

V r = I — 7^— I , X 



x-y , b > bo (M-branch) 
b 2 ~J ' * ~ I x , be(b ,b 2 ) (Bl-branch) 

(5.30) 

the critical temperature and the critical pressure are 
given by 

1 - 8xQ 2 



irv r 



2nv 2 



-in 



FIG. 18. Critical volume. The dependence of the critical 
volume v c on parameter b is depicted for Q = 1. The upper 
(black) solid line displays the M-branch. It exists for all b > bo 
and asymptotes to the RN value vrn = 2y / 6Q (upper dashed 
line). The critical point of the Bl-branch (lower solid line) 
exists for b £ (60,62), but has positive pressure only for b £ 
(60,61); the corresponding line is displayed by the red solid 
line. In this range of parameter 6 we have two critical points 
with positive pressure, one for each branch. The critical point 
of M-branch has higher pressure, temperature, and -^7^ ratio 
(see the following three figures). 



Let us note that the Bl-branch critical point has negative 
(unphysical) pressure for 



b > h = 



V3 + 273 0.4237 



6Q 



Q 



(5.32) 



Excluding this range we conclude that the Bl-branch crit- 
ical point occurs for b £ (60,61). The dependence of the 
critical quantities on the parameter 6 as well as the ratio 
p c = £^ is illustrated in Figs. [T8H2T1 The values of the 
quantities for large 6, described by the M-branch, asymp- 
tote to those of the Reissner-N6rdstrom-AdS black hole, 
given by formulas (|A14[) and (|A15[) . 

We finally mention that for 6 < 62 the black hole 
has one horizon, i.e., it belongs to the S-type. On the 
other hand for 6 > 62 (the M-branch) the black hole 
may be of S or RN type. The transition occurs when 
M(6, P c , v c , Q) = M m (b, Q). Namely, for 



6 > 63 



1.72846 



(5.33) 



the critical point occurs for the RN-type black hole, 
whereas it occurs for the S-type for 6 < 63. Obviously, it 
is the S-type black holes which possess more interesting 
properties. 

To study the possible phase transitions in the system, 
let us now turn to the expression for the Gibbs free en- 
ergy. 
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FIG. 19. Critical pressure. The dependence of the critical 
pressure P c on parameter b is depicted for Q — 1. The upper 
(black) solid line displays the M-branch, while the lower solid 
red line displays the positive pressure part of the Bl-branch. 
The dashed horizontal line stands for the asymptotic RN- 
AdS value Prn = l/(967rQ 2 ), towards which the M-branch 
asymptotes. 



FIG. 21. Critical ratio. The dependence of the critical ra- 
tio p c — on parameter 6 is depicted for Q = 1. The 
dashed horizontal is the RN-AdS asymptotic value 3/8, to- 
wards which the M-branch asymptotes. 



2. Gibbs free energy 



0.04 



0.02- 




FIG. 20. Critical temperature. The dependence of the 
critical temperature T c on parameter b is depicted for Q = 
1. The upper (black) solid line displays the M-branch, while 
the lower solid red line displays the positive pressure part 
of the Bl-branch. The dashed horizontal line stands for the 
asymptotic RN-AdS value Trn = VE / (18ttQ) , towards which 
the M-branch asymptotes. 



In order to find the Gibbs free energy of the system let 
us calculate its Euclidean action. For a fixed charge Q, 
one considers a surface integral 



1 

8tt J d- 
1 

4tt ./, 



d 3 xVhK 
d 3 xVh 



n a F ab A b 



y/1 + 



(5.34) 



The first term is the standard Gibbons-Hawking term 
while the latter term is needed to impose fixed Q as a 
boundary condition at infinity. The total action is then 
given by 



BI 



(5.35) 



where I em is given by (|5.ip , and I c represents the invari- 
ant counterterms needed to cure the infrared divergences 
[H[l3|. Thc total action reads [H (see also [H for the 
higher-dimensional case) 




(5.36) 
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FIG. 22. Gibbs free energy for b > 62. The behaviour of 
the Gibbs free energy is depicted as a function of temperature 
for fixed pressure. For b > bi the behaviour is reminiscent of 
that of the RN-AdS black hole, i.e., there is one critical point 
and the corresponding first order phase transition between 
small and large black holes for T < T c . We have set Q = 1 
and b = 1. 



The action has been calculated for fixed A and we asso- 
ciate it with the Gibbs free energy for fixed charge, 

Here, r + is understood as a function of pressure and tem- 
perature, r + = r+(P,T), via equation of state (|5.2ip . 

The behaviour of G is depicted in Figs. I221E51 It de- 
pends crucially on the value of parameter b. Namely, for 
b > 62 it is similar to that of the RN-AdS black hole, i.e., 
there is one critical point and the corresponding first or- 
der phase transition between small and large black holes 
for T < T c . The P — T diagram remains qualitatively 
the same as for the 4D RN-AdS black hole and is almost 
identical to Fig. [5] 

For b £ (61,62) one has one physical (with positive 
pressure) critical point and the corresponding first oder 
phase transition between small and large black holes. 
This phase transition occurs for T < T C1 and terminates 
at T = Tt . There also exists a certain range of temper- 
atures, T G (T t ,T 2 ), for which the global minimum of 
G is discontinuous, see Fig. [251 In this range of tem- 
peratures two separate branches of intermediate size and 
small size black holes co-exist. They are separated by a 
finite jump in G. Although there is no real phase transi- 
tion between them, we refer (for simplicity) to this phe- 
nomenon as 'zeroth-order phase transition'. For T < T t 



G 




0.02 0.03 0.04 0.05 



FIG. 23. Gibbs free energy for b £ (bi, 62). In this range of 
parameters one has one physical (with positive pressure) criti- 
cal point and the corresponding first oder phase transition be- 
tween small and large black holes, occurring for T £ (Ti , T C1 ) 
and P € (Pt,P cl ). Note that contrary to the case 6 > 62, 
small black holes (r+ — > 0) now correspond to high tempera- 
ture, T — > 00. There also exists a certain range of pressures, 
P e (P t ,Pz), and temperatures, T G (T t ,T z ), for which the 
global minimum of G is discontinuous (see Fig. I26[) and, aside 
from the first order phase transition, there are two phases of 
intermediate and small black holes separated by a finite jump 
in G. Fro T < Tt only one phase of large black holes exists, 
see P—T diagram in Fig. [2H1 We have set Q = 1 and b = 0.45. 



only one phase of large black holes exists, see P — T di- 
agram in Fig. [27] and its magnification in Fig. E51 

For 6 G (60, 61) there exist two critical points with pos- 
itive pressure. However, only the one at T = T Cl corre- 
sponds to the first order phase transition between small 
and large black holes. The other, with lower temperature, 
lower pressure, smaller radius r + , and smaller mass, does 
not globally minimize the Gibbs free energy and hence is 
unphysical. As 6 decreases the two critical points move 
closer to each other and finally merge and disappear com- 
pletely for 6 = 6 . Similar to the case 6 G (61,62), there 
is a 'zeroth order phase transition' between intermediate 
and small black holes. For T < T t only a phase of large 
black holes exists. The P — T diagram is very similar to 
the case when 6 G (61,62)- 

Finally, for 6 < 60 there is no first order phase transi- 
tion in the system and the behaviour is like that of the 
Schwarzschild-AdS black hole. 



3. Critical exponents 

In the preceding discussion, we have learned that for all 
6 > 60 there is a first order (small-large)-black hole phase 
transition in the system, which terminates at the critical 
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FIG. 24. Gibbs free energy for b G (bo,bi). In this range 
one has two critical points with positive pressure. However, 
only the one at T = T C1 corresponds to the first order phase 
transition between small and large black holes. The other 
does not globally minimize the Gibbs free energy and hence 
is unphysical. As b decreases these critical points move closer 
each other and finally merge and disappear completely for 
b — bo. Similar to the previous figure in between (Tt,T z ) there 
is a 'zeroth order phase transition' between intermediate and 
small black holes. We have set Q — 1 and b = 0.4. 



point described by (i> c , T c , P c ), given by Eqs. (|5.30l) and 
(|5.3ip . This critical point occurs for the M-branch, i.e., 
it is described by x = xi, (|5.27|) . [The other possible crit- 
ical point belonging to the Bl-branch is unphysical since 
the corresponding phase of the system does not globally 
minimize the Gibbs free energy] For large b the critical 
values expand as 



v c = 2V6Q 
V6 



216 Qb 2 v 6 4 



V6 



18Qtt 2592irQ 3 b 2 
1 



96n( 



and asymptote to those of the RN- AdS black hole (|A14I) . 
Exact values (for any b > bo) are given by (|5.30[) and 
(JOT]) and arc displayed in Figs. [T%rl2Tjl (M-branch) . The 
critical ratio p c — is displayed in Fig. [5TJ For large 
b it expands as 



1 



8 384Q 2 & 2 



(5.39) 



and asymptotes to the value 3/8, characteristic for the 
Van der Waals fluid, (|A"8|) , or the RN-AdS black hole, 



FIG. 25. Gibbs free energy for b < bo. For b < bo, similar 
to the Schwarzschild-AdS case, there is no first order phase 
transition in the system. We have set Q = 1 and b = 0.3. 




0.03 



0.038 



FIG. 26. Zeroth order phase transition. The global min- 
imum of the Gibbs free energy (displayed by the black solid 
line) is highlighted by the red thick line. We have set Q = 1, 
b = 0.45 and P = 0.46 x Pci w 1.6854 x 10" 3 G (Pt, Pz). Ob- 
viously, there is a first order phase transition between small 
and large black holes occurring at T = T\ ~ 3.3340 x 10~ 2 . 
There is also a discontinuity in the global minimum of G at 
T — To ~ 3.2799 x 10~ 2 . This separates the intermediate size 
black holes from the small black holes. We refer to this as 
a 'zeroth-order phase transition'. This type of behaviour is 
characteristic for all b G (&o,&2) and P G (Pt,Pz)- 



17 



Pel 
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0.002 - 




T 

0.06 



FIG. 27. P -T diagram for b G (61,62). The coexistence 
line of the first order phase transition between small and large 
black holes is depicted by a thick solid black line for Q = 1 
and 6 = 0.45. It initiates from the critical point (P C1 ,T C1 ) 
and terminates at (Pt,Tt). There is also a 'zeroth order phase 
transition' (dotted line) discussed in more detail in the next 
figure. 



0.0018 - 



0.0017 



0.0016- 




FIG. 28. P — T diagram magnified: zeroth order phase 
transition. The first order phase transition between small 
and large black holes is displayed by thick solid black line. 
The red solid line describes the 'coexistence line' of small 
and intermediate black holes, separated by a finite gap in G, 
indicating the zeroth order phase transition. It commences 
from (T Z ,P Z ) and terminates at (Pt,Tt). 



Let us now study the vicinity of this critical point by 
calculating its critical exponents. The entropy 5, as func- 
tion of T and V is given by 



1/3 



(5.40) 



It follows that Cy = and hence the critical exponent 
a = 0. To calculate other critical exponents we employ 
the equation of state (|5.22|) . We introduce the quantities 
p, to and t, given by (|2.25|) and (|2.28[) with q = 3, and 
series expand the law of corresponding states, which is 
again of the form (|2.27[) , to obtain 



p = 1 + At - Btui - Clu 3 + O(toj 2 ,uj 4 ) 



where 



A = 



D 



1 

3/3 c 



(5.41) 



(5.42) 



and C, given by a relatively complicated expression, 
whose exact form is not really important. For simplicity, 
we give the large b expansion 



C-— 7 of — 

~ 8 1 29 1 66 2 Q 2 + 



(5.43) 



which obviously asymptotes to the RN-AdS value 4/81. 
Following the discussion in Sec. II. C, it is now obvious 
that the critical exponents a, (3, 7, and S take the same 
values as in the (higher-dimensional) Maxwell case, 



0, p 



7 = 1, 5 = 3. 



(5.44) 



So we have shown that even for the Born-Infeld black 
holes, obeying highly nonlinear electromagnetic field 
equations, we obtain the same critical exponents as in 
the linear Maxwell case. 

Let us finally mention that the derived critical ex- 
ponents (|5.44[) differ from those of recent papers [30I — 
|33j. The authors therein study the critical behaviour of 
the Born-Infcld-AdS black hole system in non-extended 
phase space and obtain various different from us critical 
exponents which do not agree with the mean field theory 
prediction. The reason for this is that they do not really 
study the vicinity of the critical point, but rather the be- 
haviour close to the two points characterized by horizon 
radii r\ and r-i where the specific heat at constant charge 
Cq diverges. Such points, however, correspond to the 
critical point only in the limit when these radii coincide, 
i.e., r\ = 7'2 = r c , in which case the analysis in these 
papers is no longer valid. In all other cases the two stud- 
ied points do not correspond to the global minimum of 
the Gibbs free energy and hence are unphysical. Hence 
the analysis in these papers is not correct. We expect 
that a proper analysis, e.g., similar to [H| . would lead to 
the same value of the critical exponents as derived in our 
paper. 
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FIG. 29. Vacuum polarization B. The behaviour of 6 as a 
function of b and r+ is displayed. Apparently, as we approach 
the linear Maxwell case (b increases) or we are far from the 
source (r+ increases) B rapidly decreases. For r+ = and 
non-trivial b we have a finite value of B, as further confirmed 
by analysis and the following figure. 



E. Vacuum polarization 

Let us finally briefly return to the new quantity B, 

(5.45) 

which was recently computed as a conjugate quantity to 
the Born-Infeld parameter b (29[. An interesting ques- 
tion is a physical interpretation of this quantity. Since 
b has units of the electric intensity and M units of en- 
ergy, we can see from formula (|5.14[) that B has units of 
electric polarization (or more precisely a polarization per 
unit volume). Hence we refer to B 'Born-Infeld vacuum 
polarization'. This quantity is necessarily non-vanishing 
in Born-Infeld space times and is essential for consistency 
of the Smarr relation ([5.141) . 

The behaviour of B as a function of r + and b is depicted 
in Figs. Inland I5IH Apparently, as we approach the linear 
Maxwell case (b increases) or we are further away from 
the source (r+ increases) B rapidly decreases to zero. On 
the hand, for r + — > and a non-trivial value of &, the 
quantity B reaches a finite positive value, see Fig. [30] 

Using Eq. (|5.2T|) we can rewrite r+ as function of 
P, T, b and Q, thereby obtaining from (|5.45[) a 'new' equa- 
tion of state 

B = B(b,Q,P,T) (5.46) 

in the extended phase space (including b and B variables) . 
Similarly the expression (|5.37[) is now understood as G = 



B 

\ 

1- \ 

\ \ 




1 2 3 



FIG. 30. Vacuum polarization B for fixed b. The be- 
haviour of B as a function of r+ is displayed for 6=1, 0.45, 0.4 
and 0.3 (from bottom up). Note that for r+ = 0, B reaches 
the finite value. 



1.5-1 



1- 



0.5- 




FIG. 31. B — b diagram. Isobars in (Z3,6)-plane are depicted 
for various choices of P and fixed T = 3Trjv = V6/(6irQ). 
All isobars emerge from the point [6o,Bo] which corresponds 
to small black holes, r+ — > 0. The subsequent behaviour as r+ 
increases depends on the value of P. Namely, the displayed 
curves correspond to P/Prn = 500, 20, 4 and 0.2, with Prn = 
1/(96ttQ 2 ). We have set Q = 1. 



G(b, Q, P, T). It then makes sense to construct a '(B—b)- 
diagram', drawing the isobars (isotherms) for fixed Q and 
T (P) . We display this in Fig. [31] which demonstrates 
interesting behaviour that is further confirmed by the G— 
b diagram displayed in Fig. [32] We leave the analysis of 
these diagrams (and a possible critical behaviour therein) 
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FIG. 32. Gibbs free energy as function of b. G is depicted 
as function of b for various values of P and fixed T = 3Trn • 
Similar to the previous figure, all curves start from the point 
[&o,Go] and the subsequent behaviour depends on the ratio 
P/P c - We have set Q = 1 and the values of P as in previous 
figure. 



as possible future directions. 



VI. SUMMARY 

We have investigated the thermodynamic behaviour of 
a broad range of charged and rotating black holes in the 
context of an extended thermodynamic phase space. For 
Einstein-Maxwell theory this phase space includes the 
conjugate pressure/volume quantities, where the former 
is proportional to the cosmological constant and the lat- 
ter is proportional to a term that would be the Newtonian 
geometric volume of the black hole. These relationships 
hold in any dimension, and are necessary for the Smarr 
relation to be valid. We find in all dimensions d > 4 
that an analogy with the Van der Walls liquid-gas sys- 
tem holds, provided we identify the specific volume in 
the Van der Waals equation with the radius of the event 
horizon (and not the thermodynamic volume) multiplied 
by an appropriate power of the Planck length. However 
only in d — 4 does the critical ratio p c = ^r 2 -, (]2.23l) . 
equal its corresponding Van der Waals value of 3/8. The 
critical exponents also coincide with those of the Van der 
Waals system, due to a universal behaviour, Eq. (|2.27|) . 
in which the pressure is equal to a term proportional to 
the temperature divided by the specific volume plus a 
function f(v) only of the specific volume. This situation 
is a consequence of the form of the metric function in Eq. 
(|2.2|) and holds for all the black holes we consider in this 
paper, including neutral slowly rotating black holes in 
d = 4. The only exceptions are charged or rotating black 



holes in d = 3, for which no critical behaviour exists due 
to the particular form of f(v). 

We also examined the thermodynamic behaviour of 
Born-Infcld black holes. Here we found that the phase 
space needs to be further extended, to include not only 
the pressure/ volume terms from the Einstein-Maxwell 
case but also the conjugate pair (B,b). Inclusion of this 
latter pair is necessary to obtain consistency of the Smarr 
relation. The (recently computed p9j ) former quantity, 
B, we refer to as the Born-Infeld polarization because of 
its dimensionality. While we have analyzed to some ex- 
tent the behaviour of this quantity in Figs. |3"T1 and I3"2l it 
is clear that its full physical meaning remains to be un- 
derstood. We find for Born-Infeld black holes two inter- 
esting types — the RN-type and the S-type — in which the 
horizon structure is quite distinct. While for sufficiently 
large b the behaviour of the Gibbs free energy is the same 
as for the Einstein-Maxwell case, once b < b-2 = 1/(2Q) 
a second critical point emerges, and the thermodynamic 
behaviour markedly changes (though the critical expo- 
nents do not change). The most intriguing feature we 
observe is that the Gibbs free energy becomes discontin- 
uous in a certain temperature range, indicating a new 
kind of phase transition between small and intermediate 
sized black holes. We call this a zeroth-order phase tran- 
sition, in keeping with the nomenclature that a first-order 
phase transition is one in which the Gibbs-free energy 
is continuous but not diffcrcntiable whereas for a sec- 
ond order transition both the Gibbs free energy and its 
first derivatives are continuous. While the physics of this 
discontinuity remains to be understood, we expect that 
the intermediate/small transition will result in a burst 
of radiation from the black hole as a consequence of the 
sudden drop in the Gibbs free energy. 

None of the examples we considered are outside of the 
the mean field theory prediction for critical exponents. 
Evidently the inclusion of non-linear electrodynamic ef- 
fects are not sufficient to go beyond mean field theory, 
and so it will be necessary to include quantum gravi- 
tational corrections to do so. A recent attempt along 
these lines (35j posits an ansatz for the quantum cor- 
rected black hole temperature for Reissner-Nordstrom 
black holes and then proceeds to investigate how the criti- 
cal exponents will correspondingly be modified. However 
the critical exponents are inappropriately identified, and 
so the question remains open. We expect to return to 
this and other interesting subjects in the future. 
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Appendix A: Critical exponents 



The critical ratio Using these quantities and defining 



Critical exponents describe the behaviour of physical 
quantities near the critical point. It is believed that they 
are universal, i.e., they do not depend on the details of 
the physical system, though they may depend on the di- 
mension of the system and/or the range of interactions. 
In D > 4 dimensions they can be calculated by using the 
mean field (Landau's) theory, e.g., [36j . 

Let us consider a fluid characterized by the critical 
point which occurs at critical temperature T c , critical 
specific volume v c , and critical pressure P c . We further 
define 



T -T c 

t = = T - 1 



v — v c 



u-1. (Al) 



Then the critical exponents a, /3,7, <5 are defined as fol- 
lows: 

• Exponent a governs the behaviour of the specific 
heat at constant volume, 



C — T— 



« \t\~ a 



(A2) 



• Exponent (3 describes the behaviour of the order 
parameter r\ = v g — vi (the difference of the volume 
of the gas v g phase and the volume of the liquid 
phase vi) on the given isotherm 



r\ = v g - vi oc \tf 



(A3) 



• Exponent 7 determines the behaviour of the 
isothermal compressibility kt. 



kt 



1 dv 
vdP 



«|t|- 



(A4) 



• Exponent S governs the following behaviour on the 
critical isotherm T = T r : 



IP -P c | (X \V-V C 



(A5) 



Iii particular, for the fluid described by the Van der 
Waals equation in any dimension d > 3 (see appendix B) 



(A6) 



P +Z1 )(v-b)=T, 



p v T 

P =P C > V = ^ T =T C > 

we can rewrite the Van der Waals equation in the univer- 
sal (independent of a and b) law of corresponding states 



8r = (3v - 1) [p + — 



(A10) 



The critical exponents are then calculated to be (see, e.g., 



a = 0, 13 = -, 7=1, 5 = 3, 



(All) 



which are the values predicted by the mean field theory. 

Similarly, for a spherical charged AdS black hole in 4D 
we have the following equation of state [3]: 

P_I » + 2* (AU) 

V ZlTV z 7TV 

Here, P stands for the pressure associated with the cos- 
mological constant, P = — J^A, T is the the black hole 
temperature, Q its charge, and v is the corresponding 
'fluid specific volume' which can be associated with the 
horizon radius r + as 

v = 2l 2 P r+ . (AI3) 
This equation admits a critical point which occurs at 



T r = -^L , v c = 2V6Q , P 



18ttQ 



96ttQ 2 



(A14) 



Similar to the Van der Waals equation, these quantities 
satisfy 

v „, q 

(A15) 



independently of the charge of the black hole. Using fur- 
ther definitions (|A9j) . we get the universal (independent 
of Q) 'law of corresponding states' 



8t = ^(p + ^)-^- 

The critical exponents are calculated to be 



(A16) 



where we have set the Boltzmann constant k = 1, and 
a > 0, b > are constants describing the size of molecules 
and their interactions. The critical point occurs at 



a = 0, 0=-, 7=1, 6 = 3, (A17) 
and coincide with those of the Van der Waals fluid. 



8a a 
Tc= 276' Vc = 3b > Pc= 27P 

giving the universal critical ratio 

P c v c 3 



Pc 



T c 



(A7) 



(A8) 



Appendix B: Van der Waals equation in higher 
dimensions: heuristic derivation 



Let us repeat the statistical mechanics derivation of 
an equation of state of a weakly interacting gas. In d 
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number of spacetime dimensions, the partition function 
of TV interacting particles with positions and momenta 
ri , pi , i = 1, . . . , TV is given by 



dri . . . dr^dpi . . . dp Are 



-0H 



2m ' 

i<3 



(Bl) 



with /3 = kT, and Wij being the two-particle interac- 
tions. Integrating over momenta, while using the formula 



dp cxp(— P§—) = J^jr, an d expanding the inter- 
acting part of the exponential {/3 » J2 Wij), we get 



,n „ iV(d-l) 
' 2lTm\ 2 



/3 J 
'27rm\ 



f dr x ... dr N (l-/3j2w, 



TV(TV - 1) 
2V 2 



dr 1 dr 2 Wr2j , 
(B2) 



with V being the volume of the gas. Going to the centre 
of mass system, we can re-write the last integral as 



W r2 dr 1 dr 2 = —2Va , a 



1 



W(\r\)dr. (B3) 



Using large TV approximation we finally have 



/3TV 2 a^ 



(B4) 



The free energy is then defined as F = —kTlgZ, which 
taking a small, expansion writes as 



F 



- — -fcTTVlg(27TTOfcT)-fcTTVlgU- — + 0(a 2 ) 



V 



(B5) 



This thermodynamic potential obeys dF = —PdV 

d 

dv 



SdT + [idN . The equation of state then reads P = — §£, 



which in terms of the specific volume v = V/N reads 



P+ — )v = kT. 

v z , 



(B6) 



One can also calculate the entropy from S = — f^j or 
calculate the specific heats for constant volume or pres- 
sure 



C — T— 

Cv ~ 1 dT 

C — T— 



fc(d- 1)TV 



k(d+l)N t 2PN 2 
P = 2 + W « + 0(a)-(B7) 



Hence we find 

_ C p _ d + 1 



4P 



d-1 k 2 T 2 (d-l 



■a + 0(a 2 ), (B8) 



generalizing the Myer's law, n = 5/3 + 0(a) to higher 
dimensions. As a final step we mention that it is quite 



natural to generalize (|B6I) by taking into account the fi- 
nite volume of gas molecules b and hence 'deriving' the 
Van dcr Waals equation in higher dimensions: 



{ p +£) (v-b) = kT. 



(B9) 



Obviously, such an equation of state inherits all the prop- 
erties of the standard Van der Waals equation. 



Appendix C: Hypergeometric identities 

Hypergeometric function p F q is defined for a\, . . . , a p , 
bi, . . . , bp, z € C as follows 



p Fq(ai,...,a p ;b 1 ,...,b q ;z) = ^ 



3=0 



(ai)j . . ■ (a p )j z j 
(Mj ■ ■ ■ (b q )j j\ ' 



(CI) 

where the Pochhammer symbol is defined as (/x)o = 
1 and (n)j = n(fi, + 1) . . . (fj, + j - 1) for j = 1,2,... 

In what follows we concentrate on 2 Fi needed in the 
main text, for which the series converges for \z\ < 1 (and 
\z\ = 1 if 3?(c — a — b) > 0). We have the following 
asymptotic expansion as z — > 0: 

2 fJ-,---. z) =i_JL + l z 2 , 3) (C2) 

V 4 2' 4' J 10 24 v ; ' v ; 

A known integral representation for 3?(c) > 0, 3?(6) > 0, 
(and any z modulo a cut along the real axis from 1 to 
infinity), reads 

, . x r ( c ) f 1 , t b - 1 (i-t) c - b - 1 

2 F 1 (a,b;c;z)= ) ' / dt i L- , 

T(b)T(c - b) J (1 - zt) a 

which for a = 1/2 and c = b + 1 reduces to 

T 



F 1 [-,b;b + l;-z) = b dt 



Using this formula we have 



d.r 



2 -Pi 



y/T+zi 



115 4 



Integrating by parts we get 



(C3) 
(C4) 

(C5) 



dx(\Jx 4 + r%- x 2 ) 



which is used in the main text. We also find 
d.r 



1 t( U 2 



o v/x 4 + r-4 4r o 0F V4 



(C7) 
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Hence, by expanding around r = we get 



Finally, using the formula for the derivative of 2-F1 



dx 



dx 



dx 



= -^ 7= r(i) 2 -4 + 0(r 2 ). (C8) 



Integrating further by parts we get the following expan- 
sion around r — 0, used in the main text: 



° / 

(y x 4 + rj] — x )dx 



3 3 
r2 _/l\ 2 



f/.r 



2= r(i)*-r 2 r + O(r 2 ). (09) 



-^- 2 Fi(a,&;c;z) = — 2 Fi(a + 1, b + 1; c + 1; z) , (CIO) 

az c 

the recurrence relation 



2 -Fi(a+l,6; c;z) - 2 Fi(a,b; c;z) 
bz 

= — 2Fi(a+ 1,6 + l,c+ l;z) 
c 

and the fact that 

2^1(5/4, 1/2; 5/4; -z) 



(Cll) 



(C12) 



we find the following two important relations used in the 
main text: 



d A 1 5 
■ 2-f*i I -; -; -z 



1 



Vl+i 



115 

2-ri I -; -; —z 



(C13) 
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1 p ^ 1 5 
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